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Abstract 

In this paper, we derive explicit sharp two-sided estimates for the Dirichlet heat kernels, 
in C 1,1 open sets D in R d , of a large class of subordinate Brownian motions with Gaussian 
components. When D is bounded, our sharp two-sided Dirichlet heat kernel estimates hold 
for all t > 0. Integrating the heat kernel estimates with respect to the time variable t, we 
obtain sharp two-sided estimates for the Green functions, in bounded C 1,1 open sets, of such 
subordinate Brownian motions with Gaussian components. 
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1 Introduction 

It is well-known that, for a second order elliptic differential operator C on M rf satisfying some 
natural conditions, there is a diffusion process X on M. d with C as its infinitesimal generator. The 
fundamental solution p(t, x, y) of dtu = Cu (also called the heat kernel of C) is the transition density 
of X. Such relationship is also true for a large class of Markov processes with discontinuous sample 
paths, which constitute an important family of stochastic processes in probability theory that have 
been widely used in various applications. Thus obtaining sharp two-sided estimates for p(t, x, y) is 
a fundamental problem in both analysis and probability theory. 

Two-sided heat kernel estimates for diffusions on M. d have a long history and many beautiful 
results have been established. See [201 [22] an d the references therein. But, due to the complication 
near the boundary, two-sided estimates for the transition density (equivalently, the Dirichlet heat 
kernels) of killed Brownian motion in a connected open set D have been established only recently. 
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See [2TJ[22j[23] for upper bound estimates and J3TJ for the lower bound estimates of the Dirichlet heat 
kernels in bounded C 1 ' 1 connected open sets. For discontinuous processes (or, non-local operators), 
the study of their global heat kernel estimates started quite recently. See EJ [HI [151 02] an d 
the references therein. See also [3] for a recent survey on this topic. The study of sharp two-sided 
Dirichlet heat kernel estimates for discontinuous processes is even more recent. In 0, we obtained 
sharp two-sided estimates for the Dirichlet heat kernel of the fractional Laplacian A Q / 2 in any C 1,1 
open set D with zero exterior condition on D° (or equivalently, the transition density function of 
the killed a-stable process in D). 

In the last few years, the approach developed in [7J has served as a road map for establishing 
sharp two-sided Dirichlet heat kernel estimates for other purely discontinuous processes in open 
subsets of R d . In E], the ideas of [7J were adapted and further developed to establish sharp 
two-sided Dirichlet heat kernel estimates of censored stable-like processes, mixed stable processes 
and relativistic stable processes in C 1,1 open subsets of R . In [2j, a Varopoulos type factorization 
estimate in terms of surviving probabilities was obtained for the transition densities of symmetric 
stable processes in re-fat open sets. Very recently, we obtained, in [12], a Varopoulos type factor- 
ization estimate for the Dirichlet heat kernels in non-smooth open sets for a large class of purely 
discontinuous subordinate Brownian motions. We have also obtained in |12j explicit sharp two- 
sided Dirichlet heat kernel estimates for a large class of subordinate Brownian motions in C ' open 
sets. 

Things become more complicated when one deals with Levy processes having both Gaussian 
and jump components. In [ID] , sharped two-sided heat kernel estimates in C ,x open sets are 
established for Levy processes that can be written as the independent sum of a Brownian motion 
and a symmetric a-stable process. A key ingredient is the boundary Harnack principle for A + A°/ 2 
in C 1 ' open sets with explicit boundary decay rates, obtained in [13] . 

The purpose of this paper is to establish sharp two-sided Dirichlet heat kernel estimates, in 
C ' open sets, for a large class of subordinate Brownian motions with Gaussian components. 
Throughout this paper, we will always assume that S = (St : t > 0) is a complete subordinator 
with a positive drift and, without loss of generality, we shall assume that the drift of S is equal to 
1. That is, the Laplace exponent of S is a complete Bernstein function which can be written as 

0(A):=A + V(A) with V(A) := / (1 - e- xt )fi(dt), (1.1) 

where /i is a measure on (0, oo) satisfying Jq°°(1 A t)fi(dt) < oo. /i is called the Levy measure of 
the subordinator S (or of (/>). We will exclude the trivial case of St = t, that is, the case of ip = 0. 
By the definition of complete Bernstein functions, the Levy measure [i has a complete monotone 
density. By a slight abuse of notation we will denote the density by n(t). For basic facts on complete 
Bernstein functions, we refer the reader to [30]. In this paper, we will assume the following growth 
condition on /i(t) near zero: For any K > 0, there exists c = c(K) > 1 such that 

H(r) < cn(2r), re(0,K). (1.2) 

Suppose that B = (Bt : t > 0) is a Brownian motion in M. d with infinitesimal generator A and 
independent of S. Then the process X = (Xt : t > 0) defined by Xt = B$ t is called a subordinate 
Brownian motion. X can be written as the independent sum of a Brownian motion and a purely 
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discontinuous subordinate Brownian motion. The infinitesimal generator of X is 

C x ■= -(f>(-A) = A - ip(-A), 
and A) is a non-local operator. The Levy density J of X is given by 



J(x)=j(\x\) 



(4vrt) 



-d/2-\x\ 2 /it 



(1.3) 



The function J{x) determines a Levy system for X, which describes the jumps of the process X: 
for any non-negative measurable function / on R + x M. d x M. d with /(s, y, y) = for all y G R d , any 
stopping time T (with respect to the filtration of X) and any x G M. d , 



/( S ,X S _,X S 



E 7 



/ (/ f(s,X s ,y)J(X s 

JO \JR d 



y)dy ds 



:i-4) 



(see, for example, p~U Proof of Lemma 4.7] and [T5j Appendix A]). 

The function j is obviously a decreasing function on (0,oo). Using this and the fact that J is 
a Levy density (and so J R d(l A |x| 2 )j(|;c|)Gte < oo), we can easily get that, for any K > 0, there 
exists c = c(K) > such that 



j(r) < cr 



-d— 2 



for r £ (0,K]. 



;i.5) 



In fact, we have for s G (0, if], 



1 



d + 2 



d+2 



r d+1 j(s)(ir < / r d+1 j(r)(ir < 
o Jo Jo 



K 



r d+l j(r)dr —. c(K) < oo, 



from which f)1.5|) follows immediately. 

The subordinate Brownian motion X has a transition density p(£, x, y) with respect to the 
Lebesgue measure. Observe that pit, x, y) is given by p(t, x, y) = p(t, \x — y\) where 



p(t,r) 



(4vrs)" d/2 e-5iP(5 4 G ds) Vt>0,r>0. 



(1.6) 



Clearly r — > p(t, r) is monotonically deceasing. For any open set D C M. d , we will use X D to denote 
the part process of X killed upon leaving D. The process X D has a transition density pn(t,x,y) 
with respect to the Lebesgue measure on D. The density po(t,x,y) is the fundamental solution of 
C x in L> with zero exterior condition, which is also called the Dirichlet heat kernel of C x in D. 

The goal of this paper is to establish explicit sharp two-sided estimates for pi)(t,x,y) in C ' 
open sets D under the above assumptions. Throughout the remainder of this paper, we assume 
that d > 1. The Euclidean distance between x and y will be denoted as \x — y\. We will use B(x, r) 
to denote the open ball centered at x G M. d with radius r > 0. For a, b G R, a A 6 := min{a, 6} 
and a V b := max{a,6}. We also define J(x,y) = J(y — x). For any Borel set A C we will 
use diam(j4) to denote its diameter and |^4| to denote its Lebesgue measure. For any two positive 
functions / and g, f x g means that there is a positive constant c > 1 so that c" 1 g < f < eg on 
their common domain of definition. 
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We will show in this paper (see Remark 12.71 below) that under the above assumptions, for every 
T > 0, there exists a constant cq = cq(T, d, tp) > so that 

p(t,r) > c (V d / 2 e- r ' 2 / 4 * + r d l 2 A (tj(r))) 

for all (t,r) G (0, T] x [0, oo). To get an explicit Dirichlet heat kernel upper bound estimates, we 
will need to assume the following upper bound condition on p(t, r) for r <diam(D): For any T > 0, 
there exist Cj > 1, j = 1,2,3, such that for all (t,r) G (0,T] x [0, diam(D)], 

p{t, r) < Ci (t-m e -r 2 IC2t + t -d/2 A (tj( r /C 3 ))) . (1.7) 

In [16], a DeGiorgi-Nash-Moser-Aronson type theory has been established for a large class of 
symmetric Markov processes on R rf with infinitesimal generators of the form 



1,3=1 x J 



Cu{x) = i-[ aii(a?) r sr z ) + lim / ( n (y) - u ( x ))] jjkrr i\ dy ^ ( Xi 



where {<Hj{x))i<i,j<d is a measurable d x d matrix- valued measurable function on Mr that is uni- 
formly elliptic and bounded, c(x, y) is a measurable symmetric kernel that is bounded between two 
positive constants, and 3?(r) is a positive increasing function in r G (0, oo). If satisfies suitable 
growth conditions near zero and infinity, sharp two-sided estimates on the transition density of this 
class of Markov processes have been obtained in [16]. In this case, the transition density p(t,x,y) 
of such a process admits the following estimates: for any T > 0, there exist Cj > 1 such that for all 
(t,x,y) G (0,T] x R d x R d , 

p(t,x,y) > ^(r^e-^-^ + t-^ACtjCcslx-yl))) 
P(t,x,y) < ci^-^e-l^/^ + t-^A^dx-yl/ca))) 

with j(r) = r~~ d (<I>(r)) _1 . These estimates can be regarded as the counterpart of Aronson's esti- 
mates for non-local operators. When (a^) is a constant matrix and c(x,y) = c{\x — y\) is a function 
of \x — y\, the Markov process X with generator C of (II. 8p is a rotationally symmetric Levy process 
on M d with Levy measure j{\z\)dz = uj^j^un in this case, the Levy exponent of X is 

*(0 = E a *Mj + / (i - cos(^ • z)) i^nn « G Rd (i-^) 

For the subordinate Brownian motion X considered in this paper, its Levy exponent is </>(|£| 2 ), 
where cj) is defined in (jl.ip , which admits an expression of the form 

m\ 2 ) = i£i 2 + / (i - cos(e • ^)) j{ Z )dz = iei 2 + voci 2 ), (1.10) 



where J is defined in (jl.3p . When the complete Bernstein function ^ satisfies the following condition 
near infinity: there exist constants 6^ G (0, 1), a& > 0, k = 1, 2, and i?i > such that 



A" 51 ^) < ip(Xr) < a 2 \ 52 ip(r) for A > 1 and r > R h (1.11) 
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then (see |28l Lemma 3.2]) 

J( x ) x i ,.,,1 , x for \x\ < 1, 

with a positive increasing function satisfying the growth conditions in |16j for r < 1. Then one 
can use the heat kernel estimates in [16] and an argument similar to the proof of [5, Theorem 2.4] 
to show that the estimate (|1.7p holds for X for any bounded open set D with C3 = 1. 

If, in addition to (jl.lip . ip also satisfies the following condition near zero: there exist constants 
5k G (0, 1), a,k > 0, h = 3, 4, and R 2 > such that 

a 3 A fi3 V>(r) < V(Ar) < a A \ 54 i[i(r) for A < 1 and r < R 2 , (1.12) 

then (see |28|. Theorem 3.4]) 

x 1 i,^,, ix for x 7^ 
v ; |x| d $(|x|) r 

with a positive increasing function <3? satisfying the conditions in |16j for all r > 0. So it follows 
from the heat kernel estimates in [16] that the estimate (|1.7p holds for X with D = M. d and C3 = 1. 

To state the main result of this paper, we first recall that an open set D in R d (when d > 2) 
is said to be a (uniform) C ' open set if there exist a localization radius Ro > and a constant 
Ao > such that for every z G 51?, there exist a C^-function ip = cp z : IR^ -1 — >■ M satisfying 
y?(0) = 0, V<p(0) = (0, . . . ,0), llV^Hoo < A , |V(^(j;) - Vip(z)\ < A \x - z\, and an orthonormal 
coordinate system CS Z with its origin at z such that 

B(z, R ) n D = {y = (y, y d ) in CS Z : \y\ < R , yd > <p(y)}- 

The pair (Rq,Aq) is called the characteristics of the C > open set D. Note that a C 1,1 open set 
L> with characteristics (i?0;Ao) can be unbounded and disconnected; the distance between two 
distinct components of D is at least Rq. Let 5qd(x) be the Euclidean distance between x and <9-D. 
It is well-known that any C ' open set D satisfies both the uniform interior ball condition and the 
uniform exterior ball condition: there exists ro = ro(i?o, Ao) G (0, Rq] such that for any x G D with 
5dD{x) < r o and y G M d \ D with 5Qrj(y) < ro, there are z x ,2; y G dD so that |x — z x \ = 5qd(x), 
\y - z y \ = 5 9D (y) and that B(x ,r ) C D and B(y ,r ) cM d \D for x = + r (x - z x )/\x - z x \ 
and yo = z y + ro(y — z y )/\y — z y \. By a C ,1 open set in R we mean an open set which can be 
written as the union of disjoint intervals so that the minimum of the lengths of all these intervals 
is positive and the minimum of the distances between these intervals is positive. 

For an open set Dcl d and x G D, we will use 5d{x) to denote the Euclidean distance between 
x and D c . For an open set D C M. d and Ao G [1, 00), we say the path distance in each connected 
component of D is comparable to the Euclidean distance with characteristic Ao if for every x, y in 
the same component of D there is a rectifiable curve I in D connecting x to y such that the length 
of I is no larger than Ao|a; — y\. Clearly, such a property holds for all bounded C > open sets, C 1 ' 
open sets with compact complements and connected open sets above graphs of C ' functions. 

For any open set D C M d and positive constants c\ and C2, we define 

h D>cl , C2 (t,x,y) := (lA^M) ( X A ^7t) (t-^e-^-y^ + t-^ A (tJ(c 2 x,c 2 y))) . (1.13) 
The following is the main result of this paper. 
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Theorem 1.1 Suppose that X is a subordinate Brownian motion with Levy exponent 0(|£| 2 ) with 
(j) being a complete Bernstein function satisfying (jl.ip and (|1.2p . Suppose that D is a C ' open 
subset ofM. d with characteristics (Rq,Aq). 

(i) If the path distance in each connected component of D is comparable to the Euclidean distance 

with characteristic Xq, then for every T > 0, there exist c\ = c\(Ro, Aq, Xq, T, -ip, d) > and 
C2 = C2 (Rq , Ao , Ao , d) > such that for all (t,x,y) e(0,T]xDx D, 

p D (t,x,y) > ci/t D)C2j i(t,a?,y). 

(ii) If D satisfies ()1 .7j) . £/ien /or every T > 0, i/iere exists C3 = C3(i?0; Ao, T, d, ip, C\, C3, d) > 1 
suc/i i/iai for all (t, x, y) e(0,T]xDx D, 

PD(t,x,y) < c 3 h Dt c 4 ,c 5 {t,x,y), 

where C 4 = (I6C2)" 1 and C 5 = (8 V 4C3)" 1 . 

(iii) If D is bounded, then for every T > 0, there exists C4 = C4(diam(D), Rq, Aq, T, ip, d) > such 
that for all (t, x, y) £ [T, 00) x D x D, 

p D (t,x,y) > C4,e~ Xlt 5 D (x)5 D (y), 

where — Ai < is i/ie largest eigenvalue of the generator of X D . 

(iv) If D is bounded and satisfies (jl.7p . i/ien /or every T > ; t/iere exists C5 = C5(diam(D), 
Ro,A ,T,^,d,Ci,C 2 , C 3 ) > sitc/i that for all (t,x,y) G [T, 00) x D x D, 

p D (t,x,y) < c 5 e~ Al< 5 D (x)5 D (y). 

When D = B(xo,r), it follows as a special case of pH Theorem 4.5(h)] that <£(A? )/2 < Ai < 
0(Af ), where Af is the smallest eigenvalue of —A in D. It follows from the scaling property 
of Brownian motion (or Laplacian) that \^( x °' r ) = cr~ 2 , where c = c(d) is a positive constant 
that depends only on the dimension d. When D is a bounded C 1 ' 1 open set in R with C 1,1 - 
characteristics (Ro, Aq), D contains a ball of radius ro and is contained in a ball of radius diam(-D), 
where ro = ro(i?0) Ao) is such that D satisfies the uniform interior ball condition with radius ro- By 
the domain monotonicity of the first eigenvalue Ax, one concludes from above that Ai is bounded 
between two positive constants that depend only on Rq, Ao, tp, diam(D) and d. 

Note that the explicit upper bound estimates in Theorem 1 1.1 1 are established under the assump- 
tion that the upper bound (|1.7|) for p(t, x, y) holds. If, instead of (|1.7|) . we assume that there exist 
constants 5 G (0, 1) and Cq > such that the function ip in (jl.ip has the property 

^(Ar) < C 6 Afy(r) for A > 1 and r > 1, (1.14) 

we can establish the following upper bound on the Dirichlet heat kernel in terms of p(t, x, y) and 
the boundary decay terms. 

Theorem 1.2 Suppose that X is a subordinate Brownian motion with Levy exponent <^(|£| 2 ) with 
4> being a complete Bernstein function satisfying (jl.ip . (jl.2p and (|1.14p . Suppose that D is a C 1 ' 
open subset ofM. d with characteristics (Rq,Aq). 
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(i) For every T > 0, there exists c\ = ci(Ro,Ao,T,ip,d) > such that for all t G (0,T] and all 

x,y G D, 

PD (t, x, y) < ci (l A ^M) (l A ^M) p(t, s/4, y/4). 

(ii) If D is bounded, then for every T > 0, there exists C2 = C2(diam(.D), Ro,Ao,d,T,ip) > 1 such 

that for all (t, x, y) G [T, oo) x D x D, 

c^ 1 e~ x ^ 5 D {x)5 D {y) < p D (t,x,y) < c 2 e" Al ' 5 D (x)6 D (y), (1.15) 

where — Ai < is i/ie largest eigenvalue of the generator of X D . 

By integrating the two-sided heat kernel estimates in Theorem 11.11 with respect to t, we can 
easily obtain sharp two-sided estimates on the Green function Gn(x,y) := J* Q °° pr)(t,x,y)dt. For 
this, let 

^ a ( 1 A fe^y ) when(i>3, 



9D(x,y) ■■= < 



log + ^gy^ when d = 2, (1.16) 

(fe(x)fe(y)) 1/2 A ^y whend = l. 



Corollary 1.3 Suppose that X is a subordinate Brownian motion with Levy exponent ^(|£| 2 ) with 
<f> being a complete Bernstein function satisfying (jl.ljl and (jl.2p . Suppose that D is a bounded C 1 ' 
open subset ofM. d with characteristics (Ro,Aq). 

(i) There exists c\ = ci(diam(Z)), Rq, Aq, ip, d) > such that 

G D (x,y) > c 1 g D (x,y), x,y € D. 

(ii) If D satisfies (jl.7p . t/ien i/iere exists C2 = C2(dia,m(D),Ro,Ao,vp,Ci,C2,C3,d) > suc/i that 

G D (x,y) < c 2 gD(x,y), x,y G D. 

We remark that even though D may be disconnected, in contrast with the Brownian motion case, 
the process X D is always irreducible because X D can jump from one component of D to another. 
Denote by G Q D {x,y) the Green function of Brownian motion in D. It is known (see [19]) that 
G G D {x,y) x gr)(x,y) when x and y are in the same component of D, and G° D (x,y) = otherwise. 
Thus when D is a bounded C 1 ' 1 connected open subset of M a! , the estimates in Corollary 11.31 are 
exactly the same as those for Brownian motion, while our heat kernel estimates (Theorem ll.l( i)- 
(ii)) detect a short-time and short-distance region, precisely t < \x — y\ 2 < 1 and 5d{x) A<5d(?/) > Vt, 
where the jump part is the dominant term. When 0(A) = A + A"/ 2 , Theorem 11.11 and Corollary 11.31 
in particular recover the main results of [10J. 

Throughout this paper the constants ro, Rq, Ao, Ao, and Ci, i = 1, . . . , 6, will be fixed. We use 
c 1j c 2) • • • to denote generic constants, whose exact values are not important and can change from 
one appearance to another. The labeling of the constants c%,C2, ■ ■ ■ starts anew in the statement 
of each result. We use c(a, /3, ...) to indicate a positive constant that depends on the parameters a, 
(3, ... Dependence on dimension d will not be explicitly mentioned. We will use dx to denote the 
Lebesgue measure in W 1 . 
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2 Lower bound estimate 

In this section we derive the lower bound estimate on po(t, %, y) when D is a C ' open set such that 
the path distance in each connected component of D is comparable to the Euclidean distance. As 
a consequence, we also get the lower bound estimate on p(t, x, y) in IR d . We will use some relation 
between killed subordinate Brownian motions and subordinate killed Brownian motions. In this 
paper we always assume that X is a subordinate Brownian motion with Levy exponent ^(|£| 2 ) with 
4> being a complete Bernstein function satisfying (jl.ip and (|1.2|) . 

Let St be a subordinator whose Laplace exponent tp is given by (II. lj) . Then t + St is a sub- 
ordinator which has the same law as St- Assume that St is independent of the Brownian motion 
B in M. d . Suppose that U is an open subset of M. d . We denote by B u the part process of B killed 
upon leaving U. The process {Zf : t > 0} defined by Z^ = B u ~ is called a subordinate killed 
Brownian motion in U. Let qu(t, x, y) be the transition density of Z u . Denote by ( z,u the lifetime 
of Z . Clearly, Zf = B t+ g for every t E [0, C Z ' U )- Therefore we have 

pu(t,z,w) > qu(t,z,w) for (t,z,w) G (0, oo) X U X U. (2.1) 

By [TJ Proposition III. 8], for every b > 0, there exist constants To > and c > so that 

f{St <bt)> c for t < T . 

Using the Markov property of St, we can easily deduce that for every b > and T > 0, there exists 
c = c(b, T, tp) > such that 

P(S t <bt)> c for t < T. (2.2) 
These facts (with 6 = 1) will be used in the proof of the following lemma. 

Lemma 2.1 Suppose that D is a C 1 ' 1 open set in Mr with characteristics (Rq,Aq) such that the 
path distance in each connected component of D is comparable to the Euclidean distance with char- 
acteristic \q. For every constant T > 0, there exist positive constants c\ = c\{Rq, Ao, Ao, MT, tp) 
and C2 = C2(Ro, Aq, Xq) such that for all A G (0,M], t £ (0, T] and x,y in the same connected 
component of y/\D, 

Pv - XD (Xt,x,y) > Cl (Xt)-^ (l A (l A e - C2 |^l 2 /(A*). (2 . 3 ) 

Proof. Suppose that X^^x and X~ l l 2 y are in the same component, say U, of D. Let pu(t, z, w) be 
the transition density of B u . By |18} Theorem 3.3] (see also |3H Theorem 1.2]), there exist positive 
constants c\ = ci(Rq,Aq, Xo,T) and C2 = C2(Ro, Ao, Ao) such that for any (s, z, w) £ (0, 2T] x U x U, 

Puis, z, w) > Cl (l A *°Vf\ (l A ,-^ e -c 2 |,-H 2 / S . (2 .4) 

(Although not explicitly mentioned in |18| . a careful examination of the proofs in [18] reveals 
that the constants c\ and C2 in the above lower bound estimate can be chosen to depend only on 
(Rq, Ao, Ao, T) and (Ro, Ao, Ao), respectively.) By using the scaling property of Brownian motion, 
we get that, for every A > 0, t G (0, T] and x, y in y/XU, 

p^uiXt^y) = X- d l 2 pu{t,X- l l 2 x,X- l / 2 y)- 
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Thus by dZ3l) . 

P^x,,) > c l( Ar^(lAM^V 1A M^)J e -«l»-»IV M 

= dfAt)-" 2 (l A ^=p) (l A ^=p) e-°l"-»l 2 /< A, >. (2.5) 

Now we assume A € (0, M]. Recall that is independent of B and that q^ u {t : x,y) is the 
transition density of z/* u = B"^f . Note that for every < t < T and x, y in y/\U, 

9^xu{^t,x,y) = / p ^ u (s,x,y)F(Xt + S\ t G ds). 
Jxt 

So by ([21D, (|22|) and ([23]), for every < t < T, A G (0, M] and as, y in a/A^, 

r-2\t 

/•At _ 

= / ^Vac/( a ^ + s > x > vMSxt G ds) 

JO 

> C3 A (l A (At)- d / 2 e- c ^^l 2 /( Ai )p(^ < At) 



> C4 

= c 4 I 1 A 



1 A (l A S -^f^j {\t)- d l 2 e- c ^-y\ 2 l^ 



□ 



Remark 2.2 Note that the Brownian motion B in M. d with infinitesimal generator A has transition 
density 

p(t,x,y) = (4vrt)~ d/2 e" j£ ^ L , x ,yGlR d , i>0. 

Using this instead of (|2.4p and an argument similar to (but easier than) the proof of Lemma 12.11 
with A = 1 we can get that, for any T > 0, there exists a positive constant c = c(T, such that 
for all t e (0, T] and x, y in 

p(t,x,y) > c 1 t~ d/2 e~^ JL 4*- . 

Lemma 2.3 For any positive constants R and a, there exists c = c(R, a,ip) > such that for all 
zeM. d and re (0,R], 

inf F y (T B ( z ,2r) > ar2 ) > c. 

yEB(z,r) 
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Proof. By Lemma |2. 11 we have 

inf P y (T B r z2r ) > ar 2 ) > P (rB(o,r) > ar<2 ) = / PB(o,r)(ar 2 ,0,y)dy 
Z -B{z,r) JBtO.r) 



y&B{z 

> / PB(o,r){ar 2 ,0,y)dy 

JB(p,r/2) 



> I Cl 

'B(0,r/2) 



C2 / a 

'B(0,l/2) 

c 3 (R,a,ip) > 0. 



V var 2 / V var 2 / 

- d/2 A A W°)) A A W*)^ e _ C2k|2/a& 



□ 



Recall that we assume that (jl,2p holds. On the other hand, since is a complete Bernstein 
function, it follows from [26, Lemma 2.1] that there exists c\ > 1 such that fi(t) < c\fj,(t + 1) for 
every t > 1. Thus by [251 Proposition 13.3.5] and its proof, we have that for any K > 0, there 
exists C2 = C2(K) > 1 such that 



j(r) < c 2 j(2r), Vre(0,if], 



and, there exists C3 > 1 such that 



j(r) < c 3 j(r + 1), Mr > 1. 



(2.6) 



(2.7) 



Lemma 2.4 Suppose that R > and 6 > 1 are constants. Then there exists c = c(R, b,ip) > 
such that for all r G (0, 12], t G [r 2 /6, 6r 2 ] and u,v & M. d , 

PB(u,r)uB(v,r)(t,U,v) > C {t~ d/2 A (tj(u, v))). 

Proof. Let r G (0,12], t G [r 2 /6,6r 2 ] and E = B(u,r) U B(v,r). If |u - u| < r/2, by Lemma O 
(with T = b, \/A = r and D = B(0, 1)) and ([13]) . 

p E (t,u,v) > inf PB(o,r)(*,0,z) = inf PB(o,r)(r 2 (t/r 2 ), 0, z) 

|z|<r/2 |«|<r/2 
> Cl t- d / 2 (l A (l A ^ e^/* > c 3 t- d/2 > C4(tJ(«,«) A t^ 2 ). 

If |ti — v\ > r/2, we have by the strong Markov property and the Levy system of X in (|1.4p that 

p E (t,u,v) > E u p E (t - T B{u , r/ ^,X TB{ur/&v v) : T B(Mjr/8) < t,X TB(ur/s) G B(v,r/8) 



PB(u,r/8)(s,u,w) [ / J(w,z)p E (t- s,z,v)dz \ dw \ ds 

WB(u,r/8) WB(u,r/8) / / 



> i inf J(w,z 

y w&B(u,r/8),z£B(v,r/8) J J 



P« ( r B(«,r/8) > s) / 

\yB(t;,r/8) 



s, z, v)dz ds 
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> p u(rB( u , r /8) > *) ( o/ /o inf J(w,z)) [ [ p B( r/8) (t - s, z,v)dzds 

\wEB(u,r/8),zeB(v,r/8) J J J B (v,r/8) 



Mrs^r/S) > t) ( w6BKr/8 ^ eB(Vjr/8) i(l- - -I)) l MTB^r/S) > s)d 

t (F (r mr/8) > br*)) 2 ( inf j(\w - z\) 

\w£B(u,r/&), z£B(v,r/8) 



> 

> c$t ( inf j(\w — z\) ) . 

\w<=B(u,r/8),z£B(v,r/8) ) 

In the last inequality we have used Lemma [2 .31 Note that, if w G B(u,r/8) and z G B(v,r/8), then 

r R 
\w — z\ < \u — w\ + \u — v\ + \v — z\ < \u — v\ + — < (2\u — v\) A (\u — v\ + —). 

Thus using both (|2.6p and ()2.7p we have 

PE(t, u, v) > c 6 tj(\u - v\) > c 6 {tJ(u, v) A r d ' 2 ) . 
The proof is now complete. □ 

The next lemma in particular implies that if x and y are in different components of D, the 
jumping kernel component of the heat kernel dominates the Gaussian component. 

Lemma 2.5 For any given positive constants ci,C2,R and T, there is a positive constant C3 = 
C3(R,T,ci,C2,ip) so that 

t -d/2 e - r y( Cl t) < C3 ( t -d/2 A (tj{ C2 r))) for every r > R and t G (0,T]. (2.8) 



Proof. Observe that (|2.7p implies that there exist C4 > and C5 > such that 

j(c2r) > C/±e~ c '° r for every r > l/c2- (2-9) 
For r > (l/c 2 ) V (2cic 5 T) and t G (0, T], we have r 2 /(2cit) > c 5 r and 

^d^-lg-r 2 /^!*) < t -d/2-l e -((l/ C2 )V(2c lC5 r)) 2 /(2cit) 

< SUp s -^-l e -((l/c 2 )V(2 ClC5 T))V(2 Cl s) =; Ce < ^ 
0<s<T 

So by ([231), when r > (l/c 2 ) V (2cic 5 T) and t G (0,T], we have 

t -d/2 e -r 2 /( C1 i) < C6te -rV(2c 1 t) < ^g-csr < (^/^ j ( 2 1Q ) 

When < r < (l/c 2 ) V (2 Cl c 5 T) and i G (0,T], clearly 

t - d /2 e -rV(c lt ) < t ^ SU p S -^-i e -« 2 /( C ^)^ < C7 t i( c 2 r). (2.11) 

The desired inequality (j2.8j) now follows from ()2.10p and (|2.1ip . □ 
Recall that the function hr> tCljC2 (t,x,y) is defined in (|1 . 13|) . 
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Theorem 2.6 Suppose that D is a C 1 ' 1 open set in M. d with characteristics (Rq,Aq) such that the 
path distance in each connected component of D is comparable to the Euclidean distance with char- 
acteristic Ao- For every T > 0, there exist c\ = ci(i?o, Ao, \o,T,ip) > and C2 = C2(i?o, Ao, Ao) > 
such that for all (t, x, y) G (0, T] x D x D, 

p D (t, x, y) > cth D}C2tl (t, x, y). (2.12) 

Proof. First note that the distance between two distinct connected components of D is at 
least Rq. Hence in view of Lemmas 12.11 and 12.51 we only need to show that there exists c = 
c(R , A , Ao,T,V>) > such that for all (t,x,y) G (0,T\ x D x D, 

Pn(t, x,y)>c (^M A A (*M A A {{tJ{X) y)) A r d l*). (2.13) 

Since D is a C ' open set, as mentioned earlier, it satisfies the uniform interior and uniform 
exterior ball conditions with radius ro = ro(Ro, Ao) G (0, i?o]- Set To = (ro/4) 2 . Consequently, there 
exists L = L(ro) > 1 such that, for all t G (0, To] and x,y G D, we can choose G D n £?(x, L\/t) 
and G DnB(y, Ly/i) so that -B(£* , 2\/i) and B(£y, 2yt) are subsets of the connected components 
of D that contains x and y, respectively. 

We first consider the case t G (0,T ]. Note that for u G B((*, we have 

S D (u)>y/i and \x - u\ < \x - £*| + |£* - u\ < Ly/i + \ft = (L + l)Vt. 
Thus by $T3§ (with A = 1), for t G (0, T ], 

*,«)*! > ci A l) [ ( 6 -^Al) t^e-^-^du 

> Cl A l) t^/2 e -c 2 (L + l) 2 |B( £ ^)| > ^ A A (2 M) 

Similarly, for t G (0,Tq], 



p D (t/3, y, u)du > c 3 f A X l • 



(2.15) 



By the semigroup property, for t G (0, Tq 



PD(t,x,y)> / / p D (t/3,x,u)p D (t/3,u,v)p D (t/3,v,y)dudv. (2.16) 

JB(&,Vi) JB(ni,Vt) 

We consider the cases |a; — y\ > \/t/8 and |x — y\ < Vi/8 separately. 

Case 1: Suppose |a; — y\ > Vt/8 and t G (0, To]. Note that by (|2.16p . Lemma [27^1 symmetry and 
pTiH2TT5]) . 



p D (t,x,y) 



> 



lB(a,V~t) JB(ti,V~t) 



p D (t/3, x, u)p B , * /2)uB , v ^ /2) (t/3, u, v)p D (t/3, v, y)dudv 
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>c 4 / / p D (t/3,x,u)(t~ d/2 A (tJ(u,v)))p D (t/3,v,y)dudv 

JB(p v ,Vi) JB(&,Vt) 

>c 4 I inf (t' d/2 A (tJ(u,v))) ) / / p D {t/3,x,u)p D (t/3,v,y)dudv 



C^ M „ M)) )(^ A1 )(!$> A1 ). 



>c 4 c§ inf (r^A(tJ(M,t)))) P^Al R-Al . (2.17) 

\(u,tj)6B(4i,Vi)xB(S*,v^) ' 

Since |x — y| > \/i/8, we have that for (u, v) G X -B(£y, V*), 

\u-v\<\u- £\ + \£ -x\ + \x-y\ + \y- fj| + |£* - u| 
<2(l + L) v / t+|x-2/| < (l6(l + L)|x-y|) A (2(1 + L)^ + |x - y|) , 

thus using (|2.6p and (|2.7p we have 



inf (t~ d / 2 A (U(u, «))) > c 5 (^ d/2 A (t J(x, y))). (2.18) 

( U ,t,)6B(^,v / I)xi?(5t,v / I) 

Thus combining with (|2.17p and (|2.18p . we conclude that, for \x — y\ > y/i/8, 

PD (t, x, y) > c 6 A l) (^M A ^ ((t J(x> y)) A r ^/2). (2 . 19) 

Case 2: Suppose \x - y\ < y/t/8 and i G (0,T ]. Then for (u,u) G \/t) x y/i), 

\u - v\ < 2(1 + L)y/i + \x-y\< (2(1 + L) + 8^)Vt. 
Thus by (|2Mi|) . we have for every (ti, v) G x B(^ y ,Vt), 

p D (t/3,u,v) > c 7 A l) (^L A t - rf /2 e - C8 |«-^A > Cg ^/2. 

Therefore by (l2TTi1l - (l2TT6| ). for t < T , 

P o(t, X ,y) > ^1(^ai)(^Mai)^/ 2 

> c 9 c§ (^Al)(^ (2.20) 

Combining (|2TT9|) and (I2T201) we get (|2TT3|) for t G (0, T ]. When T > T and t G (T , T], observe 
that T /3 < t - 2T /3 < T - 2T /3 < (T/T - 2/3)T , that is, t - 2T /3 is comparable to T /3 with 
some universal constants that depend only on T and Tq. Using the inequality 



p D (t,x,y)> p D (T /3,x,u)p D (t - 2T /3,u,v)p D (T /3,v,y)dudv (2.21) 

instead of (|2.16p and by considering the CclSGS | X 

y\ > y/To/8 and \x — y\ < a/To/8 separately, we 
deduce by the same argument as above that (|2.13p holds for t G [To, T] and hence for t G (0, T\. □ 
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Remark 2.7 By Lemma [2.4l we have that for every T > there is a positive constant c\ = c\ (tp, T) 
such that for all (f , x, y) G (0, T] x R d x R d , 



p(t, x, y) > P B(x ^ t)vjB{y ^ t) (t, x, y) > ci d/2 A (i J [x, y)fj 



(2.22) 



Together with Remark 12.21 (|2.22j) yields the following global lower bound on p(t,x,y): For every 
T > 0, there is a positive constant C2 = 02(1(1, T) such that for all (t,x,y) G (0,T] x R d x M d , 



3 Upper bound estimate 

In this section, we derive the upper bound estimate on p£>(t,x,y) for C 1 ' open sets satisfying the 
assumption (|1.7|) , We first record a lemma, Lemma I3.1} which serves as the starting point for the 
upper bound estimate. Applying it and using (|1.7|) for po(t, x, y) on the right hand side of (|3.2[) . we 
can get an intermediate upper bound estimate for p£>(t,x,y) that has one boundary decay factor. 
This is done in Proposition 13.21 Applying Lemma 13 . 1 1 again but now using the intermediate upper 
bound estimate for p£>(t,x,y) obtained in Proposition 13.21 on the right hand side of (|3.ip . we can 
get the desired short time sharp upper bound estimate for pu[t,x,y). This is carried out in the 
proof of Theorem ll.lf ii). Recall that X is a subordinate Brownian motion with Levy exponent 
0(|£| 2 ) with (f) being a complete Bernstein function satisfying (jl.ip and (jl.2p . 

Lemma 3.1 Suppose that U±, C/3, E are open subsets ofM. d with U\, U3 C E and dist([/i, Us) > 0. 
Let U2 := E \ (U\ U U 3 ). If x G Ui and y € U 3 , then for every t > 0, 



Proof. The proof is similar to that of [10\ Lemma 3.4]. For the reader's convenience, we spell out 
the details here. Using the strong Markov property of X, we have 



p{t,x,y) > c 2 (Vrf/2 e -k-?/l 2 /(4t) + t -d/2 A ( t J( X) y))) 



(2.23) 




(3.1) 



(3.2) 



p E {t,x,y) 



E x PE(t - TU^Xr^ ,y) \T Ul <t 

Ex Pe (t - t Vi , X TUi , y) : r Ux < t, X TUi G U 2 

+ E x Pe (t ~ ru x , X TUi , y) : t Ux < t, X TUi G U 3 =: I + II . 



Clearly 
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On the other hand, by (jl.4p and the symmetry, 

Jo yJu! \Ju 3 



< 



sup J{u,z) \ / Pa;(Ti7i > s) ( / p E (t - s,z,y)dz) ds 
\ueUi,zeu 3 J Jo \Ju 3 J 

J ^ F x ( TUl > s)F v (t e > t - s)ds. 



< sup J(u, z) 



Finally 



f F x (t Vi > s)¥ y (r E > t - s)ds < [ F x ( TUl > s)ds < t A E x [ TUl ]. 
Jo Jo 

This completes the proof of the lemma. □ 

Recall that C%, C2 and C3 are the constants in ()1.7|) . 

Proposition 3.2 Suppose that D is a C 1 ' 1 open set in M d with characteristics (i?o,Ao). Assume 
that (|1.7p holds. For every T > 0, there exists c = c(C%, C3, Ro, Aq, T, tp) > such that for all 
t £ (0, T] and all x,y £ D, 

PD(t, x, V )<c(lA ^) ^- d /2 e -|,-y|V(4C 2i ) + (r d/2 A tj( | x _ y | /(4 v 2C3))) ) 

Proof. There exists ro = ro(Ro, Aq) 6 (0, i?o] such that D satisfies the uniform interior and 
uniform exterior ball conditions with radius ro. Fix T > and t £ (0, T\. Let x,y £ D. In view 
of (jl.7p . we only need to show the theorem for 5d(x) < ro\/i/(16vT) < r o/(16), which we will 
assume throughout the remainder of this proof. Choose xq G dD such that 5d(x) = \x — xq\. Let 

Ux := B(x , r Vt/(8VT)) n £>. (3.3) 

Let n(a;o) be the unit inward normal of D at the boundary point xq. Put 

lbv -£ 

Note that 5d(xi) = . Applying the boundary Harnack principle in [27] we get 
P.^ £ D\Ui) < Cl F Xl (X TUi £ D \ U,)^ < C1 ^Y X) - 

Hence 

F x (X TUi £D\U 1 )<c 2 (lA S -^\ . (3.5) 

By [271 Lemma 4.3], 

Exfa/J < c 3 Vi 5 D (x). (3.6) 
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Thus we have by (|3.5p and (|3.6p . 

P x (td > t/2) < P* > t/2) + P x (X^ £ D\Ui 

Sp(x) 

Vt 



< n ~E X [rc/J jAlj+lPx (x TUi G £> \ C/i) < c 4 f 1 A 



(3.7) 



Now we deal with two cases separately. 
Case 1: |z - y| < (VMC^ V (r /VT)) v^- By the semigroup property, symmetry and (|1.7]h 



p D (t,x,y) = / p D (t/2,x,z)p D (t/2,z,y)dz 
Jd 



< sup p(t/2,z,w) / p D (t/2,x,z)dz 

yz,weD J JD 

< 2 1+d / 2 C 1 t- d / 2 F x (T D > t/2) 

< Ci 2 1+d / 2 C lt - d / 2 (l A , 

where in the last line (I3.7P is used. Since 

\x - y\ 2 /(4C 2 t) < (d/2) V (r 2 /(4C 2 T)) < (d/2) V (r 2 /(4T)), (3.8) 

we have 

Mi,*,y) < C4 2 1+d / 2 C ie W 2 ) v ^/( 4T ))t- d / 2 e-l^^ 2 /( 4 ^*) (l A . 

Case 2: \x-y\> (j2c^V (r /VT)) V*. Let 

J7 3 := { z e £> : | z _ x | > \ x — 3/|/2} and f7 2 := -D \ (E/i U U 3 ). (3.9) 

Since 

\z — x\ > - — > ° for z 6 I/3, 



dist(f7i, C/3) > and, for u e U\ and z £ U; 



3. 



|u — 2 1 > |je — z| — \x — Xq\ — \xq — u\ > \z — x\ — ~~/= ^ \ z ~ x l/2 ~>\x — y|/4. (3.10) 



Thus, 



sup J(u,z)<c sup j(\u - z\) < c 3 j(\x - 7/1/4). (3.11) 

uel/i,zet/ 3 («,«):[«-«[>i[x-i/| 

if ze U 2 , 

2 I y I 

— y\ > \x — y\ + \x — z\ > \z — y\ > \x — y\ — \x — z\ > — - — . (3.12) 
We remark here that up to this point, we have not used assumption (jl.7p yet in this proof. 
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Observe that, for any f3 > 0, the function f(s) := s d l 2 e is increasing on the interval 
(0,2/3/d]. By dHZJ), (|3321) and the observation that t < \x - y\ 2 /(2dC 2 ), 

sup p(s,z,y) < Ci sup ( s -d/2 e -\z-y\ 2 /C 2 s + s -d/2 AsJ ( z /c 3 ,y/C 3 )) 

s<t,zeu 2 s<t,zeu 2 

< Ci sup ( S - d / 2 e-l 2 -^ 2 / C2S + S J(z/C3,y/C 3 )) 

s<t,\z-y\>\x-y\/2 

< C lS nps- d/2 e- lx - yl2/ ^ C2s) + dtj(\x -y\/(2C 3 )) 

8<t 

< dr^e-^ 2 ^^ + c 5 (r d/2 A ijflz - y|/(2C 3 ))), (3.13) 

where in the last line (jl.5p is used. In fact, since \x — y\ > (ro/VT)\/t, by (jl.5p . 

/ f \ l+d/2 

tj(k-y|/(2C 3 )) < tj((k-y|Ar )/(2C 3 ))<c 6 — t" d / 2 

\ |x — y| z A / 

< c 6 (T/r 2 ) 1+d / 2 t- d / 2 , (3.14) 

where C6 > depends only on C 3 . 

By the same argument as that used to get (|3.5p . we can apply the boundary Harnack principle 
in [27] to get 

^ X {X G U 2 ) < c 7 ¥ xl (X € tf 2 )|^- < c 8 ^. (3.15) 
Applying <^2il, fl3~TT|) . (l3~T3]l and (^T5l) . we obtain 



p D (t,x,y) 



<c 9 (t^/ 2 e -l^l 2 /(4C 2t ) + r m A tj - ()x _ y |/(2Cs))) ^ + c 10 tj(|z - y|/4)^ 

< Cll ^-^/2 e -|-J/| 2 /(4C 2 t) + t - d /2 A _ y y {2C3 v 4)) ^ 

where in the last line (|1.5p is used (see (|3.14p ). This combined with (jl.7p completes the proof of 
this proposition. □ 



Proposition 3.3 Suppose that D is a C 1,1 open set in W d with characteristics (Rq,Aq). Assume 
that (|1.7p holds. For every T > ; there exists c = c(Ci, C 3 , Ro, Aq, T, -0) > suc/i that for all 
t e (0, T] and all x £ D 

P,(r >t)<c(lA^M). 

Proof. Fix T > 0. By Proposition 13.21 and (|2.23p we have that for every < t < T and x, z in Z), 

PD(t, x, z) <d (l A ^) ^-d/^-I^P/CiW + r d/2 a (tj (|x — z|/(4 V 2C 3 ))) 

( Sd(x)\ , 
<c 3 II A 1 p(i, c 2 x, c 2 z), 
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where c 2 := (y^ V 4 V (2C3))" 1 . Thus 

<5£)(a;) 



:(td > t) = J^p D (t,x,z)dz <c 3 ^lA ^7=^) f D , 



p(t, C2X, c%z)dz < C3 I 1 A 



□ 



Proof of Theorem ll.lt (i). This has already been established in Theorem [2? 

(ii). Fix T > 0. There exists ro = ro(Ro,Ao) £ (0, Ro] such that L> satisfies the uniform interior 
and uniform exterior ball conditions with radius ro. Let t 6 (0,T] and x,y £ D. By Proposition 
13. 2\ (jl.7p and symmetry, we only need to prove (ii) for Sd(x) V Sjj(y) < ro\/t/(16vT) < r o/(16), 
which we will assume throughout the remainder of the proof of (ii) . 

The proof of (ii) is along the line of the proof of Proposition 13.21 but using the estimate from 
Proposition 13.21 for the upper bound estimate of pu(t,x,y) on the right hand side of (|3,ip rather 
than using (|1.7j) . Define U\, xq and x\ in the same way as in the proof of Proposition 13.21 (see 
(|3.3|) - (|3.4|) ). and consider the following two cases separately. 

Case 1: \x — y\ < ( y/8(d + 1)6*2 V {tq/\/T) \ y/i. By the semigroup property, symmetry and 



Proposition 13.21 

p D (t,x,y) = I p D (t/2,x,z)p D (t/2,z,y)dz 



D 



< suppD(t/2,y,z) / p D (t/2,x,z)dz 

\z£D J JD 

< Cl ^/ 2 (lA^M)p.(r D >t/2) 



< cit 



~ d ' 2 { l A ^7r) i 1 A ^7T 



< C2 t- d / 2 e -l-fl 2 /(i6C 2 t) ( , . ) J , A ™ j cuO) 

where in the third inequality, Proposition 13.31 is used. 

Case 2: \x - y\ > (y/8(d + 1)C 2 V (r /VT)) y/i. Define U 2 and U 3 as in QSXfy . Note that (I37TTT) 



^)( 1A «$>). 



holds. Moreover, for z G C/2, as in (j3. 12[) . 

3\x-y\/2 > \z-y\ > \x-y\/2. (3.17) 

Observe that, for any j3 > 0, the function /(s) := s _ ^ +1 ^ 2 e _/3//s is increasing on the interval 
(0,2/3/(d + 1)]. By Proposition EjZ] (instead of CLZD), (13~TT)) and the observation that t < \x - 
y\ 2 /(8(d + l)C 2 ), 

sup PD (s,z,y) < c 3 sup ( s -d/2 e -l^l 2 /( 4C ^ + ^ A 5j(|z-y|/(4V2C 3 )))^^ 

s<t,z&U 2 s<t,z£U 2 V s 

< c 3 5 D (y) sup ( a -(*fi)/2 e -[*- tf [7(4C a .) + _ y|/( 4 V 2C 3 ))) 

s<t,|*-vl>[a!-l/[/2 

< c 3 fe(y)(sup S -( d+1 )/ 2 e-^-^ 2 /( 16C2S ) + Vtj(|x - y\/(8 V 4C 3 ))) 
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< C4 ^) ( t -d/2 e -|*-y| a /(l6Q.t) + ( t -d/2 A _ 3,1/(8 v 4C3))) , (3.18) 

where in the last line we used an argument similar to that in (|3.14p . On the other hand, by 
Proposition 13.31 we have 



I F x (t Vi > s)F v (t d > t - s)ds < [ P x (t d > s)F y (T D > t - s)ds < c 5 f 
Jo Jo Jo 



Sd(x) 5 D (y) 
— F as 



f 1 1 

= c 5 5 D (x)5 D (y) dr = c e 5 D (x)5 D (y). (3.19) 

Jo \Jr(l — r) 

Combining ([53]) . ([BTTT]) . ([3TT5|) . (I3TT81) and (l3TT9|) all together, we conclude that 

PD(t,x,y) 

< c 7 (t-d/2 e -\x-yf/(mc 2 t) + {t - d /2 A mx _ y | /(g v 4C , 3))) ^ S D (x)5 D (y) 

+c 8 tjfls - 2/1/4) 

< c 9 (t-^ e -N-y| 2 /(i6C 2t ) + (t -d/a A tj{]x _ yU{8 v 4C3))) ^ S D (x)S D (y) 

= c 9 (t-^ e -\x-y\ 2 /^t) + {t - d /2 A y(|x _ y|/(8 v 4C3))) ^ A A L A ; 

where in the second inequality f| 1 . 5 1) is used (see (|3.14p ). This combined with ()3.16j) and Proposition 
13.21 completes the proof of (ii). 

(hi) and (iv). We first note that the path distance condition is satisfied in any bounded C 1,1 open 
set D with Ao depending only on Rq, Aq and diam(D). Thus by (i) and (ii), it suffices to prove (iii) 
and (iv) for T = 3. 

In view of (jl.7p . the transition semigroup {Pf,t > 0} of X D consists of Hilbert-Schmidt 
operators, and hence compact operators, in L 2 (D;dx). So has discrete spectrum {e~ Afc *; k > 1} 
arranged in decreasing order and repeated according to their multiplicity. Let {(f>k,k > 1} be the 
corresponding eigenfunctions with unit L 2 -norm, which forms an orthonormal basis for L 2 (D;dx). 

Clearly, 

/ (lA6 D {x))<k(x)dx<\Dn<k\\i^m)<\D\ 1/2 - (3-20) 
Jd 

By using the eigenfunction expansion of pr>(f, x, y) = J2T=i e_Afci ' ^k{x)(j>k{y) > we get 

f (1 A6 D (x))p D (t,x,y)(lA6 D (y)) dxdy = Ve" 1 ^ ( [ (1 A S D (x)) <f> k (x)dx\ . (3.21) 

JDXD \JD J 

Noting that \k is increasing and \\fW2 = f { z )^k{z)dz) 2 , we have for all t > 0, 

tAi / /1 a x_ ('^X 2 j„. ^ — iAi 



(lAfc(x))p D (t,i,j/)(lAfo(j/)) dxdy < e" iAl / (1 A 5/)(a;)) dx < e |Z?|. (3.22) 
DxD Jd 

On the other hand, by Theorem ll.in i) and (|3.20p . there exists c\ > so that for every x £ D, 

<pi (x) = e Xl / p D (l,x,y)(f>i(y)dy 
Jd 
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< cie Al 



(lAfo(x)) I (lA5 D (y))My)dy < c^D^ 2 (lAfc(x)). (3.23) 
Jd 

It now follows from f)3.21j) that for every t > 0, 

(1 A 5 D (x))p D (t,x,y) (1 A 5 D (y)) dxdy 



'DxD 

> e- tXl (J {c^D^/^^xfdx^j =c^ 2 |D|- 1 e-(* +2 ) Al . (3.24) 



For t > 3 and x,y £ D, we have that 



p D (t,x,y) = p D (l,x,z)p D (t - 2,z,w)p D (l,w,y)dzdw. (3.25) 

JOxD 

By Theorem I l.lf ii) and (|3.22p . there exist q > 0, i = 2, 3, so that for every t > 3 and x,y £ D, 

p D {t,x,y) 



< c 2 (lA5 D (x))(lA5 D (y)) / (1 A 5D(z))pD(t - 2, z,u>) (1 A 5d(vj)) dzdw 

J DxD 

< c 2 \D\e- x ^-V (1 A fo(z)) (1 A 5 D (y)) < Cge^ 1 ' (1 A 5 D {x)) (1 A 5 D {y)) . (3.26) 

By ()3.25p . Theorem 12.61 the boundedness of D and ()3.24p we have that there exist q > 0, i = 4, 5, 
so that for every t > 3 and x,y £ D, 

PD(t,x,y) 

> c 4 (lAj(dmm(D))) 2 (lA5 D (x))(lA5 D (y)) [ (1 A 5 D (z)) p D (t - 2, z,w) (1 A 5 D (w)) dzdw 



I DxD 

> c 5 (lAj(diam(Z)))) 2 l^l" 1 (1 A 5 D (x)) (1 A 5 D {y)) e~ iAl = c 6 (1 A (1 A 5 D {y)) e~ tXl . 

The theorem is now proved. □ 

Let (£, J-) be the Dirichlet form of X on L 2 {M. d \ dx). It is known that (£, F) is a regular Dirichlet 
form on L 2 (M. d ;dx) with core C^(R d ); see [3]. Moreover, for -u G C*( 



£ (u, u) := / Vu(x) • Vv(x)dx + / (u(x) — u(y)) 2 J(x,y)dxdy (3.27) 



and J 7 := C c 1 (M d ) £l C W 1 ' 2 (R d ) = {/ € L 2 (R d ;cfo) : £(/,/) < oo}. So we have the following 
Nash's inequality: 

2 +4/ "<ci/ |Vn(x)| 2 dx- ||/||^ < c 2 £(f,f)\\f\\t /d for/eJT. ( 3 .28) 
It follows then 

x, y) < c 3 t~ d/2 for t > and x, y € M d . (3.29) 
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Proof of Theorem 11.21 (i) There exists ro = ro(i?o> Ao) € (0, Rq] so that D satisfies the uniform 
interior and uniform exterior ball conditions with radius ro. Fix T > 0. We claim that there is a 
constant c$ > so that 

p D {t,x,y) <c p(t,\x-y\/A) A lj for every (t, x, y) £ (0, T] x D x D. (3.30) 

In view of (|1.6|) . po{t,x,y) < p(t, \x — y\) < p(t, \x — y|/4). So it suffices to prove (|3,30p when 
5d{x) < rQ\/t/{l(S^/T) < ro/16. Define U\, xq and x\ in the same way as in the proof of Proposition 
[331 (see ([313]) (E3])). Hence ([33]H(l3/F]) and (|3T5j) hold. We now prove (EOUl) by considering the 
following two cases. 

Case 1: \x — y\ < (ro/v^V^- By the semigroup property, symmetry and (13.29j) . 



p D (t,x,y) = / p D (t/2,x,z)p D (t/2,z,y)dz 
Jd 

< sup p(t/2,z,w)\ / p D (t/2,x,z)dz 
yz,weD J Jd 

< c lt -^P x (r D > t/2) < c 2 t- d ^ (l A , 

where in the last line (13 .7p is used. Since \x — y\ 2 /(64i) < rg/(64T), we have by Remark E7J 



PD 



(t,x,y) < C2 e r V^ T h- d l 2 e-\ x -y\ 2 l^ (l A 5 -^\ < c 3 p(t, \x - y|/4) (l A 



8 D (x) 



Vt 

Case 2: \x - y\ > (r /VT)Vi. Define U 2 and U 3 as in <^M). Note that (IBTTT]) and (f3T2l hold. 



Observe that by (|1.6|) 

sup p(s,z,y)< sup p(s, z, y) < supp(s, \x — y\/2). (3.31) 

s<t,zeU 2 s<t,\z-y\>\x-y\/2 s<t 

By [281 Lemma 3.1], (fTTIj) implies 

j(r) < c 2 r- d ip(r- 2 ) < cic 2 ^(l)r~ d ~ 25 for all r G (0, 1]. 

Thus under assumption (|1.14p . according to [IBJ Theorem 1.3], the parabolic Harnack inequality 
holds for the subordinate Brownian motion X. Extend the definition of p(t, r) by setting p(t, r) = 
for t < and r > 0. For each fixed x,y E M. d and i > with \x — y\ > (ro/ 'y/T)y/i, one can easily 
check that (s,w) h-> p(s, |w; — y|/2) is a parabolic function in (— oo,T] x B(x, (ro/v / T)v / i/4). So 
by the parabolic Harnack inequality from |16[ Theorem 1.3], there is a constant c 3 = C3(V0 — 1 so 
that for every t £ (0, T], 



Hence we have 



supp(s, |rr - y|/2) < c 3 p(i, |x - y\/2). 

s<t 



sup p(s,z,y) < c 3 p(t, \x - y\/2). (3.32) 

S<t, Z&J2 



Applying ([32D, ([M]), (13XTT) . (I3T5|) . (13321) . we obtain 



PD(t,x,y) < f x [X T eU 2 )\ sup p(s,z,y) +1^^] sup J(u,z) 
v 7 \s<t,zel/2 / Wf/i,zeC/ 3 
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< c 5 p(t,\x-y\/2) 5 -^ + c 5 tj(\x-y\/4) 6 -^. (3.33) 



Since \x — y\ > roy/i/yT, by (jl.5p . 

/ t \ 1+d/2 
tj(\x-y\/4) < tj((\x-y\Ar )/4)<c A (- ^— ^ r d ' 2 

< c 4 (T/r 2 ) 1+d / 2 ^/ 2 . (3.34) 

Thus by the monotonicity of the transition density and (|3.33|) and (|3.34p we obtain 



PD(t,x,y) < c 5 p(t,\x-y\/4)^ + c e (t~ d / 2 Atj(\x-y\/4)) 
< C7p(t,\x-y\/4) 



< c 7 p(t,\x-y\/4)( S -^-M), 



^ »u^j ^ 

Vt 
'A 

where in second inequality (|2.22p is used. 

Combining these two cases establishes the claim (|3.30p . Thus by the semigroup property and 
the symmetry of pD(t,x,y) in x and y, we conclude from (|3.30j) that for every t G (0, T] and 
x,y € D, 



p D (t,x,y) = I p D (t/2,x,z)p D (t/2,z,y)dz 

\ A (l A J pit/2, \x - z\/4)p(t/2, \z - y\/4)dz 

\ A (l A £ p(t/2, \x - z\/4)p(t/2, \z - y\/4)dz 



< c? 



< c 2 ^ 

<- -(^^)(iA^)KM.-,l/4). 

(ii) The lower bound in (j!.15j) is Theorem ll.lf iii). The proof of the upper bound in (I1.15P is 
the same as that of Theorem ll.lf iv) , the only difference is that we use part (i) of this theorem and 
(|3.29p . instead of Theorem I l.lf ii). so that Pd(1, x, z) < c\(l A 5d(x))(1 A 5d(z)) and pz)(l, u>, y) < 
cx(l A 6 D (w))(l A 6 D (y)) . □ 



4 Green function estimates 

In this section we give the proof of Corollary 11.31 

Proof of Corollary 11.3b Put T := diam(L>) 2 . Recall that gn{x,y) is defined in ()1.16|) . By an 
argument similar to that for [7\ Corollary 1.2], one gets that (see the proofs of [291 Theorem 5.0.8] 
(for d = 1, 2) and [24, Theorem 6.2] (for d > 3) for details) 

Jo ( lA ^7r) (^^) t - d/2 ^ CllX - Vl2/tdi + f e-^ S D (x)S D (y)dt * g D (x,y). 
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Thus, since D is bounded, by Theorem 11.11 (i) and (hi), we have Gr>(x,y) > C2go{x,y), which 
proves Corollary ll.3f i). 

When the bounded C ' open set D satisfies (|1.7p . by (jl.5p and Theorem 11.11 (ii) and (iv), we 
have 

OM, ,) , . (, D (x. ,) + [ (i a (l a *f ) (r- A ^) * 

Therefore to prove Corollary II. 3f ii) it suffices to show that 

T ' 1A ^7r) ( 1a (^ A F^) * <- 



o 

By the change of variable u = \ x ~^ , we have 

cT 

y/i J V v 7 * / V I s -2/1 



o 



u^ 2 A n- 3 ) ( 1 A 2&eW) ( 1 A A^M') du . (4 . 2 ) 



k - yl ' 2 J\x- y \ 2 /T ^ ' V k-2/l / V 1^ — 2/1 

Since for every x,y £ D and r > 0, 



V \ x -y\J \ \ x -y\J F-y| 



we have 



< 







-2 




1 




X 


_ y |d- 


-2 




1 




;r 


_ y |d- 


-2 




1 




;r 




-2 



oc 



u^ 2 " 2 A n- 3 ) (l A ^ffi) (l A vg^M') du 



\x - y\ 



oo 



F-2/1/ V f-2/1/ 



oo 

„.-2 



F - 2/1/ V F - 2/1/ 

a; - 



(1) When d > 3, by (@3D, we have 



1 ^/ 2 - 2 A n- 3 l f 1 A ^iM) ( i A ^*°M ) du 



|» - y\ d 2 J\x- y \2/T ^ ' V ' k-2/l / V k-2/l 

F - 2/1 V F - 2/1/ V F ~ 2/1, 



o 



< ' f 1A fefejfeM). (4.5) 

~ d-2 |x-y| d ~ 2 V k-y| 2 1 



Combining (pE2"|) , (jO|) and (g3D , we arrive at (|4TT]) for d > 3. 
For the other cases, we dehne 



(4.3) 



< ^(IaW. (4-4) 



n := Mf^M2/). (4.6) 
F ~~ 2/1 
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Clearly l/n > \x - y| 2 /diam(L>) 2 = \x - y\ 2 /T. 
(2) Suppose d = 2. By (|4.3|) we have 
1 



\x - y\ 



d-2 



< z 1 u- 1 fiA ug f (g) y icfa 

|x-y| 2 /T V \x-y\ z J 

1 _ rl 

u ll \u>i/u }du+ / lt l{ u <l/ uo }d'!i 
\x-y\ 2 /T J\x-y\ 2 /T 



= log 

Thus by ([13]) and (^71) . 



/ 1 |x — i/| 2 

«o V 1 + u — A 1 - J" 

Vuo r 



< 



1A^) flA^) + log(. Vl) + nof^A 

f - y|/ V f - y\J \ u o 

( 1 b - yl 2 

1 A n + log(u V 1) + uq — A 1 - J" 

1 A n + log(n V 1) x log(l + n ) = log (l + ^ ,(r) ^ ( "> 



\x - 



\x - j/| 5 



This proves (j4. 1 j) for d = 2. 

(hi) Lastly we consider the case d = 1 . By (|4.3p and (|4.6 
1 



k - y| 



d-2 



|*-»| a /T 



u d/2-2 



A-" 3 ) (lA ^ ] ) (lA ^^ } ) 



|*-tf| a /T 



|x - y\< 



\x-y\[ u 3/2 l {u>1/uo} du+ / n n 1/2 l{ u <i/ U0 }dn 

\x — y| 2x ' 



= \x-y\ ^2 ((n V l) 1 / 2 - l) + 2n ^(n V l)" 1 / 2 - ( 
Thus by P~2j) . and the last display, we have 

< \x - y\ (1 A n ) + |x - y| M (n V l) 1 / 2 - l) + n ^(n V I)- 1 ' 2 - 

x |x - y| (nf A n ) = (M*)Mv)) 1/!! A ■ 
This proves (14. 1|) for d = 1. 



2 
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